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Abstract 

We evaluate the static QCD potential to two-loop order. Compared to a previous 
calculation a sizable reduction of the two-loop coefficient a<i is found. 
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1 Introduction 



The static potential of (massless) QCD has recently been calculated to two loops |]]. Being 
a fundamental quantity, it is of importance in many areas, such as NRQCD, quarkonia, 
quark mass definitions and quark production at threshold. While the one-loop contribu- 
tion and the two-loop pole terms have been known for a long time @, || |J, the two-loop 
constant a 2 (cf. Eqs. (fl)-flHD) was found only recently PJ. The fermionic parts of this 
coefficient were confirmed numerically in ||, taking the m q — » limit of a calculation 
involving massive fermion loops. The aim of this work is to evaluate a 2 analytically, using 
an independent method. 

The static potential is defined in a manifestly gauge invariant way via the vacuum 
expectation value of a Wilson loop [f|. || . 

V(r) = - lim^ — In (trP exp (ig jf dx^A^ . (1) 

Here, T is taken as a rectangular loop with time extension T and spatial extension r, and 
is the vector potential in the fundamental representation. 

In a perturbative analysis it can be shown that, at least to the order needed here, all 
contributions to Eq. ([I]) containing connections to the spatial components of the gauge 
fields Ai(r, ±T/2) vanish in the limit of large time extension T. Hence, the definition can 
be reduced to 

V peIt (r) = - Jim I ln(trTexp (- jf ^)) > ( 2 ) 

where T means time ordering and the static sources separated by the distance r = \r — r '\ 
are given by 

j;{x)=ig5 lxQ T a [8{x-r)-8{x-r')} , (3) 

where T a are the generators in the fundamental representation. In the case of QCD the 
gauge group is SU{3). The calculation will be carried out for an arbitrary compact semi- 
simple Lie group with structure constants defined by the Lie algebra [T a ,T b ] = if abc T c . 
The Casimir operators of the fundamental and adjoint representation are T a T a = Cp 
and f acd f bcd = CA$ ab - tr(T a T b ) = Tp5 ab is the trace normalization, while n/ denotes the 
number of massless quarks. 

Expanding the expression in Eq. (g) perturbatively, one encounters in addition to 
the usual Feynman rules the source-gluon vertex ig5^T a , with an additional minus sign 
for the antisource. Furthermore, the time-ordering prescription generates step functions, 
which can be viewed as source propagators, analogous to the heavy-quark effective theory 
(HQET) 
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Concerning the generation of the complete set of Feynman diagrams contributing to 
the two-loop static potential, there are some subtleties connected with the logarithm 
in the definition (0). All this is explained in detail in [M pi, so we only list the relevant 
diagrams here (see Figjl]). Note that the aforementioned papers are based on the Feynman 
gauge, while we use general covariant gauges, resulting in an enlarged set of diagrams. 



2 Method 

The method employed in this work can be briefly summarized as follows: 

• All dimensionally regularized (tensor-) integrals are reduced to pure propagator 
integrals by a generalization of Tarasov's method M. The resulting expressions are then 
mapped to a minimal set of five scalar integrals by means of recurrence relations, again 
generalizing || as well as @ to the case including static (noncovariant) propagators^]. 
These two steps are implemented into a FORM [|Il|] package. Thus, we constructed our 
method to be complementary to the calculation in M , assuring a truly independent check. 
At this stage, one obtains analytic coefficient functions (depending on the generic space- 
time dimension D as well as on the color factors and the bare coupling), multiplying each 
of the basic integrals. 

• The basic scalar integrals are then solved analytically. Expanding the result around 
D = 4 (which is done in both MAPLE |12| and Mathematica [13| considering the com- 



plexity of the expressions) and renormalizing, one obtains the final result to be compared 
with 0. 

Important checks of the calculation are the comparison of the pole terms of individual 
gluonic diagrams given in ||, the gauge independence of appropriate classes of diagrams, 
the confirmation of cancellation of infrared divergencies, and the correct renormalization 
properties. 



3 Renormalization and result 



The renormalized quantities in D = 4 — e dimensions are conventionally defined by 

V = ^V R , -^ = Z^a R , (4) 



1 This algorithm will be described in detail elsewhere fiofl . 
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where the subscript R denotes the renormalized quantities. The factor Z is assumed to 
have an expansion in the renormalized coupling, Z = 1 + clrZ '(e) + a R Z" (e) + and 
we choose to work in the MS scheme, related to the MS scheme by the scale redefinition 
^ = tfei/A-K . 

The needed counterterms read explicitly 

z' =--#> , z"_ = \pl--Px. (5) 

ms e MS e e 

Here, the coefficients of the Beta function are defined by the running coupling, fi 2 d^ a R = 

-Poa 2 R - Pxa 3 R - A) = t Ca ~ I and |8i=fCS- AC F T F n f - f C A T F n f . The 
results of our calculation yield indeed Eq. (|5|). 

As a further check on the pole terms, the vertex and gluon wave function renormal- 
ization constants (Z% and Z^ 1 , respectively) have been extracted separately from the 
diagrams. They depend on the gauge parameter £ and agree with the ones given in [14|. 



The renormalized potential now reads Vp ert (r) = / -p^k exp(iqr)V (q 2 ) , with 

v„ 2) = r m + t_c,jp\ + ..\ ( 6) 



q2 ms I ms ms I q2 J ms ms I q2 



where 



Ci ms^ = 01 +A)ln(z) , (7) 
c 2 _(x) = a 2 + (3 2 \n 2 (x) + (f3 1 + 2f3 a 1 )ln(x) (8) 



and 



oi = y^i-y 2>n/ , (9) 
/4343 , 2 7T 4 22 \ 2 /1798 , 56 



° 2 = (i62- +4 ^-T + y^J^-l^r + T c(3)) ^ 



- 16C(3)J C F T F n f + ^T 2 nj . (10) 



As it has to be, the coefficients prove to be gauge independent. Comparing our two-loop 
result for a 2 with [[[]], we find a discrepancy of 2tt 2 in the pure Yang-Mills term (oc C\). 
This amounts to a 30% decrease of a 2 for the case of n/ = 0, and a 50% decrease for rij = 5 
(for 577(3)), which is the case needed for ti threshold investigations. This difference can 
be traced back to a specific set of diagrams, as outlined below. 

The origin of the discrepancy is Eq. (14) in the second paper of [0J. To explain the 
crucial point, let us introduce some notations first: We have two types of denominators, 



Di = D(ki) = -k , stemming from gluon, ghost and fermion propagators, and Si = S(ki) = 

i 

- fc ^ +jg , with v = (1, 0), stemming from the source propagators. The loop momenta are 
k\, k 2 , k% = k\ — q, k^ = k 2 — q, h = h — k2, where q = (0, q) is the external momentum. 
We abbreviate the integration measure as J { = // / ^H> > while products of propagators 
will be written like DiD 2 = Di 2 etc. 

Adding the diagrams in question gives (neglecting the color factors) 

2£ + 5K + X = li (^235^1125 + ZWll25 + D 235 S ± 



122; 



Jl I2 ^ 1 



/ 145'S'll25 + -02355*1225) 

= jJ 2 D U5 S U2 (S 5 + S- 5 ) , (11) 

where the identity 5*1 5*2 = S 5 (S 2 — Si) (compare |], §4) was used for the last term of the 
first line, and the trivial exchange of loop variables ki «-> k 2 was done in the last term of 
line two. In the last line, Sg = S(—k 5 ). One then obtains 

S 5 + S~ 5 = } : , , + ' 



(ho - k 20 ) + ie -(ho - k 20 ) + is 

- - {km -Z> + e> - -*"<*»-*»)■ < 12 > 
Hence, contrary to the assumption in QXJ] , the sum of the integrands in Eq. (|TTD reduces 
to a delta distribution multiplying the remaining propagators. 

Now, considering the color traces as well as the gluon-source couplings, one gets as a 
contribution to the bare static potential (for simplicity, we use the Feynman gauge here 
to make the point clear) 



diag.a6 + diag.b3 = - 9 - C F C\ ( ]g + - 9 - C F C\ ]g 

,6 „6 



-^<*X-^ca(X + Z + X 



While in []l] the latter term was discarded, we evaluate it in D = 4 — e dimensions to give 



(D-5) 

' i + 0(c). (13) 



32tt 2 q 

Note that the factor of (£) — 4) in the numerator cancels the single pole in the scalar 
integral, such that only the constant part of a 2 is affected by this discussion, while the 
pole terms are not changed by the omission of this term. Hence, dividing out the overall 
factor (~ ( 16 ^2^ 6 ^ ); we identify the 2t[ 2 C\ difference with respect to ]1[ Q. 
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q^ 

We thank M. Peter for checking this result. 



The static potential can be used for a definition of an effective charge, which is 
conventionally called a v . Defining V(q 2 ) = —C F 16ir 2 a v /q 2 , one can use the knowledge of 
the three-loop coefficient /3 2 MS |T3[] to derive the corresponding coefficient in the V-scheme 
from Eq. (§). While /3 and (3\ are universal, one finds 

PI = PF - axPi + (a 2 - a?)A) ( 14 ) 
44vr 2 llvr 4 242 




C(3) C\ 



12 9 

- y + ^ C(3)) C\T F n f + Q + 2 -f ((3) ) C A T F nj 



, : , C(3) J C A C F T F n f + 2C F T F n f + - - C (3)j C F Tj^ (15) 

The new value for a 2 leads, for SU(3) and = 5, to a 50% decrease of /3J compared to 
the formula given in [p]]. 

Summarizing, we have re-calculated the two-loop static potential by a method com- 
plementary to the approach in |l[]. We have developed an algorithm which enables us to 
work in general covariant gauges throughout. Confirming the fermionic contributions to 
the two-loop coefficient a 2 , we find a substantial deviation in the pure gluonic part of a 2 . 
The source of the discrepancy could be identified. 

We would like to thank W. Buchmuller, M. Spira, T. Teubner and M. Peter for 
valuable discussions and correspondence, respectively. 
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y^y ylPZx £\ y^N, STV^ 




b3 cl c2 c3 c4 c5 c6 dl 




d2 el e2 e3 fl £2 f3 




g2 hi hlg hlf h2 h3 il i2 




il.bg il.bf il.c il.d il.e il.f 



Figure 1: Classes of two-loop diagrams contributing to the static potential. Double, 
wiggly, dotted and solid lines denote source, gluon, ghost and (light) fermion propagators, 
respectively. A blob on a gluon line stands for one-loop self-energy corrections. 
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